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Abstract 

We formulate a generalization of the Glauber formalism for hadron-nucleus and nucleus- 
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small distances contribute to the Pomeron structure; and (ii) the triple Pomeron vertex 
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as the survival probability of large rapidity gap processes and inclusive observables, both 
for hadron - nucleus and nucleus- nucleus collisions. Our approach suggests saturation of the 
density of the produced hadrons in nucleus-nucleus collisions, the value of the saturation 
density turns out to be large and depends on the number of nucleons in the lightest nucleus. 
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1 Introduction 



Although QCD is beheved to be the microscopic theory of strong interactions it has led to little 
progress in the theoretical understanding of "soft" interactions at high energies. At present we 
only have a phenomenological description of such processes where the main postulate of this 
description is the Pomeron - a Reggeon with a trajectory intercept ap(0) close to unity, 

ap{t) = ap{0) + a'p{0)\t\ ; ap(0) = 1 + A, ; A<1. (1.1) 

We have two comments on the Pomeron hypothesis: 



1. The first one is negative. There is a paucity of ideas and/or examples of how a Reggeon 
such as the Pomeron can be justified theoretically. As far as we know, the only theoretical 
model, where the Pomeron appears naturally, is a 2 + 1 dimensional QCD which can 
hardly be considered as a good approximation to reality. The entity closest to the Pomeron, 
which follows from perturbative QCD (pQCD) is the BFKL Pomeron 0. In the following, 
when dealing with "soft" interactions, we shall also refer to the BFKL Pomeron, but only 
with regard to the general structure of pQCD approach at high energy. 



2. The second comment is positive. In spite of all theoretical uncertainties the Pomeron has 
been at the heart of high energy phenomenology over the past three decades, providing 
a good description of all available experimental data 0,0, p. Consequently, we cannot 
ignore the Pomeron hypothesis as well as the numerous attempts to find a selfconsistent 
approach based on this postulate. 



Even if we accept the Pomeron hypothesis, there still remains a lot of hard work to find the high 
energy scattering amplitude that includes the Pomeron exchanges, as well as the interactions 
between 

Pomerons. The goal of this paper is to find a solution to this problem for hadron-nucleus and 
nucleus-nucleus high energy coUisions. Our approach is based on two main ideas: (i) a new 
sufficiently hard scale for the Pomeron structure and (ii) a specific parameter (ka) for the Pomeron 
interaction with a nucleus suggested by Schwimmer many years ago, 

= A-'-^ ■ U(^f - 1}. (1.2) 

itRa a So 

where Qpn is the Pomeron-nucleon vertex, G^p is the triple Pomeron vertex and A is defined in 
Eq. ( p..l[ ). A and Ra are the atomic number and the radius of a nucleus, which are defined by 

Ra = A^- 17GeV-'^. (1.3) 
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1.1 A new scale for the Pomeron structure 



Our key assumption is a new scale for the Pomeron structure, i.e. a sufficiently large mean trans- 
verse momentum of partons in the Pomeron { < pt >> 1 GeV ). We list, first, the experimental 
and phenomenological observations that support our point of view (see also Ref. 



1. Experimental elasting scattering data can be fitted using a very small value for the slope of 
the Pomeron trajectory ( see Refs. §, [§, §, § ): = 0.25GeV-^ < a'^ = IGel/-^, 
here is the slope of the Reggeon trajectory ( see Eq. ( p..l|) ); 

2. The experimental ^ slope of single diffraction dissociation, with final state secondary 
hadron system with a large mass, is approximately two times smaller than the slope for 
the elastic scattering. The description of the diffractive dissociation processes is closely 
connected to the triple Pomeron vertex, G^p ( see Fig. |l| ). A small slope leads to a small 
proper size of this vertex. As a first approximation we can assign a zero slope to the 
triple Pomeron vertex, as this provides a reasonable description of the experimental single 
diffractive dissociation data. 




3. The HERA data |10| on diffractive photoproduction of J/\E' and deep inelastic scattering ( 
DIS), show that the t- slope for elastic diffractive dissociation ( 7* + p J/\E' + p ) is 
larger than the t-slope for the inelastic channel ( 7* + p J/\E' + A^*, where A^* is the 
nucleon exitation ). These data provide direct experimental support for the idea that there 
are two typical scales for the size inside a hadron, namely, the distance between quarks and 
the proper size of the constituent quarks in the Additive Quark Model( AQM )[|rT]. In the 
AQM the size of the constituent quark is closely related to the typical transverse momentum 
of partons in the Pomeron. 
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4. The main idea of the AQM that gluons inside a hadron are confined in a volume of a smaller 
radius ( Rq ~ 0.1/m <C -R^ ~ Ifm ), is a working hypothesis which is included in the 
standard Pomeron phenomenology that describes the experimental data @]. 



5. The HERA experimental data ^ have confirmed the theoretical expectation [12| that the 
density of partons in a hadron at low x reaches a high value. In high density QCD we 

This saturation scale means that the 
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expect that a new saturation scale appears 
average transverse momentum of a parton becomes large in the region of high density and, 
therefore, it confirms our assumption regarding a new scale associated with the Pomeron. 



The consequence of assuming a new hard scale in the Pomeron leads to: 



• The Pomeron exchange as well as all vertices for the Pomeron-Pomeron interaction, can 
be viewed as delta functions in impact parameter ( ht ) representation. Indeed, based on 
our assumption, we can neglect both ap(0) and the proper size of all Pomeron-Pomeron 
vertices, and consider the Green's function for the Pomeron exchange, as well as the ver- 
tices, as constants in momentum representation, or in other words, as delta functions in 
6t-represent at ion : 



Gspih) = G3P(0)5(2)(6^) ; G,p{bt) = G,p{0)6^'\bt) . (1.5) 

Consequently, in our approach, all the bt dependence is concentrated in the Pomeron - 

hadron vertices which depend on the size of a hadron. We use the simplest Gaussian 
parametrization for the Pomeron-Nucleus vertex, i.e.: 



exp (1.6) 



where the notation is the same as in Eq. (|r|). Combining Eq. (lOD with Eq. (|L|) one can 
obtain that for nucleus (Ai) - nucleus (^2) scattering, the single Pomeron exchange has the 
form ( at fixed bt ): 

J dH[ J dH\gpAAb't) P{s,b't-V\) gpAAbt-V\) = (1.7) 



D2 D2 I ^ K exp --2- • — ■ exp . 

nRi^ ■ nRi^ J \ RaJ ^ Sq J \ R^^ 

glj,-A^-A2 fs\^ ( hi 

' — ■ exp ' 



n{Rl^ +R%) \soJ ^""V R\, + R%, 



4 



We can use pQCD to estimate the values for the vertices of the Pomeron-Pomeron interac- 
tions, since the typical scale ( the Pomeron radius ) is small. Simply counting the number 
of as in the pQCD diagrams gives ( see Refs. § § H H ) 

gPN oc as ; A oc asN^ ; (1.8) 
G-sp oc a|A^c ; G^p cx a| . (1.9) 

To select the vertices of the Pomeron-Pomeron interaction we use the main principle of 
the large ( number of colours ) approximation, that has been formulated by Veneziano 



et al. |T^. We sum separately in each topological configuration the leading diagrams, 
considering Ncas ~ 1, while as -C 1. Applying these rules to Eq. (|1.8| ) and Eq. ( |1.9|) 
we see that we can take into account only the triple Pomeron interaction and neglect more 
complicated vertices. 



1.2 The order parameter «:a 

The problem of summing all triple Pomeron interactions looks hopeless without an additional 
small parameter that aids us in classifying the numerous Reggeon diagrams, which describe 
Pomeron interactions ( see Ref, []TB| for details ). Such a parameter is ka given by Eq. ( |1.2D 



and suggested by Schwimmer|0] more than thirty years ago, 

i^A = A ■ ■ -r{{ — ) - 1} • 1-10 

To clarify the meaning of ka we calculate the cross section for diffractive dissociation of an 
incoming hadron in the region of large mass for the hadron-nucleus collision. This process is 
illustrated in Fig. ^ where the incoming hadron is shown at the bottom of the diagram. Using 
Eq. ( p..4|) and Eq. ( |1.5|) , the cross section can be written as: 



a^^(h + A ^ M + A) = I dy J h J dH'^glM) 9PN{{bt - b',)') G,p e^^^^-y) + ^y 



{^Riy-''^" Jo 

AgpN AY Agp^Gsp 1 _ 1 



AgpN AY 

ttR^ 



e"^^ ■ = gPN9PA{ht = Q)P{Y) ■ KA , (1.11) 



^A 



where Eq. ( |1.6| ) was used and we introduce the notation: Y = /n(s/so), y = ln{M / sq) and 
P (Y) = e^^. From Eq. ( |1 . 1 1| ) we see that the parameter ka now indicates how large or small 



the cross section of diffractive dissociation is, in comparison with the total cross section given by 



Pomeron exchange ( see Eq. (1/7) ). The principle idea of the selection rules for Reggeon diagrams 
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is to sum all diagrams of the order k^, neglecting the diagrams that are proportional to gj^j^ or 
k'^. In other words, we consider the following set of parameters: 

ftA ~ 1 ; ftjv < 1 ; ■ /tjv < 1 ■ (1-12) 

9pn 

Below we show how this set of parameters helps, and the type of selection rules it leads to, both 
for hadron-nucleus and nucleus-nucleus interaction. Eq. ( |1.12D only holds for the case of heavy 
nuclei with A 3> 1. We consider the experimental fact that G^p/{gpN- ^) is rather small {k. 1/4 
- 1/8), which follows from measurements of aei{pp) and a^^{pp), as supportive for our approach 
( see also section 2.6 for details ). 



1.3 The Glauber-Gribov Approach 

As we have mentioned, the main goal of this paper is to find a natural generalization of Glauber- 
Gribov approach ||17|| ||18|| to hadron - nucleus and nucleus-nucleus interaction at high energies. In 



this subsection we recall the Glauber-Gribov approach for the nucleus - nucleus interaction (see 
Fig. 2 ). Fig. 2-a shows that in the Glauber-Gribov approach we view a nucleus as a collection of 



h h h h hhhhhhhhhhhh 

^» • • — A A A — — A m a — ^ ■ * — 




Fig. 2-a Fig. 2-b 

Figure 2: Nucleus - Nucleus scattering: The Glauber-Gribov approach (Fig. 2-a) and the first corrections 
(Fig. 2-b) to the Glauber-Gribov approach. Figures are taken fron Ref. |J^/ . 

almost free nucleons which could interact elastically and/or inelastically with a nucleon in another 
nucleus. In this approach the 

unitarity constraint 

2Imaei{s,bt) = \aei{s,bt)\^ + Gin{s,ht) , (1.13) 

has a simple solution 



aei{s,bt) = 




(1.14) 



Gin{s,bt) = 1 - exp{-Q{s,bt)) , (1.15) 
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assuming that the elastic amphtude af,i{s,bt) is pure imaginary at high energy. In Eq. ( 1.13 ) 
Gin{s,bt) denotes the contribution of all inelastic processes. The advantages of the Glauber- 
Gribov approach is that we can express the cross section of a nucleon-nucleon interaction in terms 
of an arbitrary real function Q{s, bt), which is known as opacity. Indeed, we can use Eq. ( |1.7| ) for 
rewriting the opacity in in the form: 



n 



A1-A2 



S{\ht 



b'. 



Sib't) d%d%' 



;i.i6) 



where cr^^^{s) is the total cross section of the nucleon-nucleon interaction and S{bt) is the 
distribution of nucleons in a nucleus S{bt) = J dzp{z, bt) where p{z, bt) is the nucleon density in a 
nucleus. In Eq. ( [L.7D we used a Gaussian parametrization for S{bt) to simplify our calculations. 
Formulae given by Eq. ( 1.14 ) and Eq. ( 1.15|) with Q defined in Eq. ( 1.16|) have been studied in 
detail two decades ago P, |T|, |2|, |l|, | 



Pq , |29| , |30| , |3l| , |32[| and they provide 
the first natural estimates for all observables in hadron-nucleus and nucleus - nucleus collisions. 



4 25, 2e, 



Using the AGK cutting rules ^ 
simplify our calculation. 



we can often reduce the number of Reggeon diagrams and 
A well known example of such a simplification is the inclusive cross 
section which can be described by a single diagram of Fig. 3 instead of the whole series shown in 
Fig.2. 




Figure 3: Single inclusive cross section for Nucleus - Nucleus scattering in the Glauber- Gribov approach. 
The corresponding cross-section is 

^A^^ ^ A,.A..'-^^^. (1.17) 

dy dy 

Eq. ( |1 . 1 ?[ ) is our point of reference for all further calculations. We will show that the interaction 
between Pomerons, that have been neglected in the Glauber-Gribov approach, change the result 
of Eq. ( |1.17D drastically reducing the A-dependence to 

M^LZM = 4 X 4 X . (1.18) 

dy dy 
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Fig.2-a shows that in the Glauber-Gribov approach we have neglected the intermediate processes 
of diffraction dissociation ( see Fig.2-b). In this paper we will take them into account replacing 
diffractive processes by large mass diffraction described by the triple Pomeron vertex ( see Fig.4 

)■ 



1.4 Motivation and structure of the paper 

Our primary goal in this paper is to solve the problem of high energy behaviour within the 
Pomeron framework. We are hopeful that our phenomenological approach will provide informa- 
tion as to what one can expect to learn from high energy soft interactions. We will concern 
ourselves here only with the problem of hadron-nucleus and nucleus-nucleus scattering. In these 
cases we have an additional simplification as we are dealing with systems with a high density of 
partons. We trust that matching our approach to that of QCD for high parton density systems 



12| will provide a bridge between the microscopic theory and the physics of large cross sections. 

Our second motivation is a practical one. RHIC is due to start operating shortly and we believe 
that it will be much easier to understand new phenomena such as quark-gluon plasma production, 
or the observation of the saturation scale if a reliable phenomenological model of the "soft' 
nucleus-nucleus interaction is available. Our approach is just a natural generalization of the 



Glauber-Gribov [0 approach for the nucleus-nucleus interaction, where we take into account 
not only the rescattering of the fastest partons, as was done in the Glauber-Gribov approach, but 
also the interaction of all partons with the target and the projectile. 

The structure of the paper is the following: In the next section we consider the hadron-nucleus 
interaction. The high energy amplitude for this reaction has been calculated by Schwimmer ^ 
but we recall the result of this calculation which we will need for the nucleus-nucleus collision. 
We then calculate the diffraction dissociation cross section as well as the survival probability for 
the hadron-nucleus reaction. The main properties of the inclusive production are discussed. At 
the end of the section we present a phenomenological application to exemplify the values that 
we are dealing with. The content of this section is not completely new and many properties of 
hadron-nucleus interaction have been studied for decades ( see Refs. |19|, 0, ^ 



26], ^ |28|, |2^, ^ ). However, the equation for the diffraction dissociation processes as well 
as the calculation of the survival probability of the large rapidity gap processes are described for 
the first time. 

The third section is the key section of this paper. In it we develop our technique for dealing 
with nucleus-nucleus scattering at high energy. We present the analytic formulas, as well as 
phenomenological estimates, for the total ion-ion cross sections at high energy, their elastic and 
diffractive dissociation cross sections and an estimate of the survival probability for large rapidity 
gap processes. In section four our results are summarized and discussed in light of the new 
experiments at RHIC. 



2 The Hadron - nucleus interaction 
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2.1 Selection rules for Reggeon diagrams 



In this subsection we demonstrate that Eq. ( |1.12|) leads to a selection rule: the main contribution 
to the hadron- nucleus amplitude given by the "fan" diagrams of Fig. ^ with the incoming hadron 
in the handle of this fan. 




Figure 4: "Fan" Pomeron diagrams. 



As has been shown in Eqs.( |TTT^), ( |1 . 14|) and ( |1.15|) we can introduce an arbitrary real function for 
the opacity Q{s, bt). This function is very convenient since we can use it to rewrite all observables, 
i.e.: 



(7 el 



1 



exp 



Q{s,bt 



exp 



Q{s,bt 



da, 



el 



dt 

f{s,t = ~q') 



I dHt{l - exp{-n{s,bt))) ; 
7r|/(s,t)|2; atot = 47r/(s,0); 
^ [ d%t e^^-^* a,iis,bt) . 

ZTT J 



(2.1) 

(2.2) 
(2.3) 
(2.4) 
(2.5) 



We now calculate the contributions of the Reggeon diagrams to the opacity Q. Neglecting the 
triple Pomeron vertex G^p we have only a single Pomeron exchange given by Eq. ( p..7| ) with Ai = 1. 



n{YM) 



9PN ' ^ AY 



7rR\ 



(2.6) 



where we take the integral over b[ in Eq. (|1.7| ) and neglect the value oi Rn as it is much smaller 
then Ra. Substituting Eq. ( p.6| ) into Eq. ( 1.14 ) and Eq. ( 1.15|) we obtain the Glauber formula 
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for the hadron-nucleus interaction [17, 18| • At first order with respect to G^p we have two "fan" 
diagrams of the type of Fig. |I] with diffraction of the incoming hadron or nucleus. It is easy to 
evaluate both these diagrams using Eq. ( |1 . 1 1|) . It turns out that their contributions to fl are 



n{h + A M + A) 
n{h + A h + M) 



gpN{bt - h[) gpA{bt)e^^ KA{ht) ; 
gpN{ht - h[) gpA{bt)e^^ KN{ht - h[] 



(2.7) 



where KA,N{bt 



KA,Nexp{- 



Accordingly to the rules of Eq. ( p.. 121) we neglect VL{h + A 



h + M) and, therefore, at this order we have only a "fan" diagram with two Pomerons attached 
to the nucleus. We have more diagrams of order G\p ( see Fig.^ where both up and bottom lines 
could be incoming nucleus ). 




Figure 5: "Fan" diagrams of order G\p with and without loops. 

The "fan" diagram, in which three Pomerons attach to the nucleus ( the second diagram of Fig.|^ 
) has a contribution to which is equal to 

^fan{Glp) = gPN gPAibt)e^'^ K\{bt) ; (2.8) 

while all other diagrams are much smaller since: 

1. The "fan" diagrams, where three Pomerons are attached to the hadron, are proportional to 
gpN gpAipt)e^^ n\{ht — b[) and, therefore, this diagram gives a much smaller contribution 
than Eq. (|0| ). 
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2. The diagrams, in which two Pomerons are attached to the hadron and two Pomerons are 
attached to the nucleus, also have a tree-like structure. However, its contribution is of the 
order of Qpn gpAipt)^^^ i^n i^A which is smaller than the contribution of Eq. (|2.8| ). 

3. The loop diagram ( the first in Fig.^] ) is proportional to Qpn QpAipt)^^^ f^lj ^"^^ '^^^ ^e 
neglected. 

Repeating this procedure for higher order of G^p we will obtain the selection rules that have been 
mentioned above, namely, only "fan" diagrams of Fig. ^ contribute to the value of the opacity 

n. 

2.2 The Hadron- nucleus amplitude at high energy 

To calculate the hadron-nucleus amplitude we have to sum the "fan" diagrams of Fig. ^ 




Figure 6: Hadron- Nucleus amplitude at high energy. 

It is easy to sum "fan" diagrams using the equation shown in Fig. |^. 
This equation has a simple analytic form: 

S {Y, h) = gpA {h) e^^ - G,p ■ ^ e^^^"^) ■ (y, h) dy , 

JO 

where VL{Y,ht) = gp^ S{Y,ht). Taking the derivative with respect to Y and using Eq. 
obtain a simple differential equation 



dS{yAh') 



A-S{y,bt)-Gsp ■ S^Y,b 



dY 

with initial conditions which follows directly from Eq. (|2.9| ): 

S{Y = 0,bt) = gpAibt) 



t) , 



(2.9) 



we 



(2.10) 



(2.11) 
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Figure 7: The graphic equation for the hadron-nucleus amplitude. 



Rewriting Eq. ( |2.10[ ) in the form (7 = Gsp/A): 



A-dY 

we obtain the solution to Eq. ( |2.12| ) and Eq. ( |2.11 



'^^^^'^':^ = S{yM)-7-SHY,W), (2.12) 



c(Yh]- (^""^ 9PA{ht) _ e'^^gpA {h) 

Finally, 

n cf\^h\ 9pn gpAih) /o1/l^ 

nfan{Y,bt) = gpNS{Y,bt) = (vh \ ^ ■ ^2-^^) 

Ha {j- jOt) -\- i 

2.3 Single diffraction dissociation 

In this subsection we derive the equation for the diffractive dissociation cross section of the 
incoming hadron. The general formula for the diffractive cross section can be obtained directly 
from Eq. ( |1.14| ) and Eq. ( |1.15| ) and it is 

a^^{h + A^M + A)= J d%te-^'-'''"^ D{Y,y,bt) , (2.15) 

where D{Y,y,bt) is the cross section for diffraction dissociation of the "fan" diagram at fixed 
bt- Here, Y = ln{s/so) and Y — y = ln{M'^/so). For the function D{Y,y,bt) we can write the 
equation which is shown in Fig.^ 

This equation generates all diagrams that contribute to the function D, but before doing so we 
comment on each term of the equation separately: 
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p 




Figure 8: Pictorial form of the equation for function D(Y,y,bt). 

1. We want to calculate the function D {Y, y, bt) on the l.h.s. of the equation. This function is 
a cross-section for the single diffractive dissociation generated by "fan" diagrams , with the 
gap in rapidity from to y and production of particles in the rapidity interval (Y, y) with 
Y -y = InM'^. 

2. The first term on the r.h.s. of the equation, is our simplest diagram of Fig.|l] where we 
substitute the full amplitude S{y,bt) for the "fan" diagrams (see Eq. ( p.9|) and Fig.|^ ) 
instead of the two Pomerons which are attached to the nucleus. The coefficient 2 comes 



from the AGK cutting rules |32[ since in the upper part of the diagram we have the cut 
Pomeron. In our equation this term plays the role of the initial condition and we obtain 
the full series of the diagrams for the diffractive cross section by iterating this term in the 
equation. 

3. The second term on the r.h.s. is written with the help of our function, D{y',y,bt). The 
coefficient 4 arises due to (i) we have a cut Pomeron in the upper part of diagram and (ii) 
there are two possible types of this diagram, which we can get by interchanging 5* {y, bt) 
and D{Y,y,bt). 

4. The third term on the r.h.s. describes the possibility of having a diffractive production 
from more than one Pomeron. The coefficient 2 from the cut Pomeron is cancelled by the 
coefficient 1/2, due to the symmetric form of the diagram ( we have {y',y,bt) in it ) 
giving an overall coefficient of 1. 



In writing the equation we have used the AGK cutting rules. We follow Refs. [^] and which 
establish the relation between different production processes including diffractive dissociation. 
It is easy to see that the iteration of the equation of Fig.|| yields a series for the diagram for 
diffractive production shown in Fig.|^. 
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Figure 9: The first terms obtained from the iteration of the equation of Fig. 8. 

This equation can be written in the form: 

D{Y,y,bt) = 2-gp^-GsP-e''^^-^^-Sfiy,bt) 

4 ■ gp_, ■ GsP r dy' e^^y-y'^ ■ S (y', h) ■ D (y', y, h) 

-'y 

+ gPN-G,p r dy't^'^^'-y'^ -Diy'.yM)- T dfD{y\y^^) . (2.16) 
Jy Jy 

We introduce a new function, SD (F, 6^), which is the cross section of the single diffraction of 
the incoming hadron in the hadron system with mass (M) smaller than ln(M^/so) Y — y 

SD{Y,yM) = r dy'S{Y,y'A) , (2.17) 
-^y 

with the initial condition 

SD{Y,Y,bt) = . (2.18) 
Integrating Eq. ( |2.16| ) we obtain an equation for SD (Y, y, ht) 

SD (r, y, h) = 2- gpN ■ G^p ■ T dy' S^ {y\ h) ■ e^^^'^^') 



^■gPN- G,p r dy' e^'-^-y'^ ■ S {y', h) ■ SD {y' , y, h, 
Jy 
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+ gpN ■ Gsp r d y' e^(^-^') • SD {y', y, h 



(2.19) 



The derivative of Eq. ( p^.l9 ) with respect to Y gives 
dSD {Y,y,bt) 



dAY 



SD{Y,y,bt) +gpM-l-{2S{Y,bt) - SD{Y,yM 



- 2 ■gpN-l-S^Y,bt) 

where 7 = G^p/A. 

Subtracting Eq. (|2.12| ) for 5* in Eq. ( p.20|) we obtain 

dF{Y,y,h) 



(2.20) 



A-dY 

where the function F is defined as 



F{Y,y,b,) - ^F'{Y,y,k 



•t) , 



F{Y,y,bt) = gpMT{Y,y,bt) = Qpn ■ {2S{Y,yM 



SD {Y,yM)} 



(2.21) 



(2.22) 



We obtain for F the same equation as for S, however there is a difference in the initial and 
boundary conditions, which are 



F{Y,YM) 
FiY,0) 



AY 



2gpNS{Y,bt 

o gpNgpAjbt) e 

2ka{YM) + 1 



gpNgpAiPt) e 

ka{YM) + 1 



Ay 



(2.23) 
(2.24) 



To understand Eq. ( p.24|) we need to recall that F{Y, 0) denotes the total inelastic cross section 
F(y, 0, 6t) = (jj„(F, 6t) due to the unitarity constraint of Eq. ( [1.13| ). On the other hand we 
can calculate the inelastic cross section using the AGK cutting rules which lead to ain{Y,bt) = 
2S{X,2ka). 

The solution of Eq. ( p.22|) with the initial and boundary conditions of Eq. ( |2.23| ) and Eq. ( p.24| ) 
is 

^ gpNgpA{bt)e^^ 



^2KA{Y,b,)-KAiy,h) + 1 • 
Using Eq. ( p. 221) we can find the function SD (Y, y, bt) defined in Eq. (|2.17| ) 

gpNe^^ ■ {KA{Y,bt) - HA{y,bt)} 



SD {Y,y,bt 



{2KAiY,bt) - nAiy,bt)} {KAiY,bt 



1} • 



(2.25) 



(2.26) 



^From Eq. (|2.25|) and Eq. ( p.l5|) we can obtain the total ( integrated over mass ) cross section for 
single diffraction production: 

CT^^ (Y) = (2.27) 
J rf^fo^e-^/"-™ 5D(F,0,6t) = J dHt{2gp^SiY,h) - F (F, 0, 6^ } e^^^""™- 



15 



Using Eq. ( ^125 



) we can also calculate D {Y, y, ht) in Eq. (|2.15 



), which is equal to 



|2KA(F,6t) - KA{y,bt) + Ir 



Eq. ( p.27| ) and Eq. ( p.28| ) solve the problem of single diffractive cross section in our approach. 



2.4 Survival Probability of Large Rapidity Gaps 

A large rapidity gap (LRG) process, is any process in which no particles are produced in sufficiently 
large rapidity region |Q. The simplest example of such process is the production of two large 
transverse momentum jets with LRG between them {Ay = — ^2! ^ 1 ) in back-to-back 
kinematics ( pa ~ —pt2 )'■ 

p + p — > (2.29) 

Mi{hadrons + jeti ((?/i,pti) } + [LRG{/\y)\ + M2{hadrons + jet2 {{yi,pti) }■ 

We believe that the exchange of colourless "hard" Pomeron is the only way to describe this 
production making it a unique source for experimental information regarding pQCD at high 
energies. Assuming factorization the cross section for this reaction can be written as 

ajet = f {Ay,yi + y2,Pti,Pt2) = -^i^'' (x^,Pti) ■ F^"^^ (x^,Pt2) ■ o^hard (ptux^x^s) , (2.30) 



where Fj;^'^ is the probability to find a parton with = ^^e^' in the proton and ahard denotes the 
cross-section of the "hard"parton-parton interaction at sufficiently high energies. This "hard" 
process is due to the exchange of a "hard" Pomeron. We need to multiply Eq. (|2.30|) by the 
damping factor < > to obtain the correct answer for the LRG cross section. This factor 
< \S\'^ >, gives a probability that no partons with x > from one proton, and no partons with 
X < x"^ from the other proton will interact with each other inelastically. Therefore Eq. ( p.30| ) 
should be rewritten in the following form: 

/ {Ay, y = yi + y2,PtuPt2) =< \S\' > •/ {Eq. (2.30)) . (2.31) 

This damping factor < \S\^ >, which is called survival probability of the LRG gap processes [04| , 



has been calculated in eikonal-type models [BO] for hadron - hadron collisions. We proceed to 



calculate < jS'p > for a hadron - nucleus collision in our approach, which allows us to go beyond 
the eikonal - type models. 

The general formula for the cross section of LRG processes is 



^LRG 



J dP ht e-^/-(^'^*) L (F, yi, ^2, ht) , (2.32) 
where L {Y,yi,y2,ht) is the cross section of LRG processes induced by "fan" diagrams. We need 



to calculate the "fan" diagrams of Fig. 10 to estimate L {Y, ht). 
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Figure 10: The Survival probability diagram. 



We now consider Fig. |T^ . The upper part of this diagram, from rapidity y2 to rapidity F, is our 
usual function Gpjq ■ S{Y — 1/2) but without gpAipt). We introduce the function R{Y — 1/2, z) = 
S{Y — ?/2, gpA{bt) —>■ z) , which is a generating function for the number of Pomerons at rapidity 
y2 in the diagrams of FigjTO|. This means that the functions C{Y — y2) in series 

R{Y-y2,z) = '£Cn{Y-y2)z^ (2.33) 

can be interpreted as the probabihty to have n Pomerons at rapidity ?/2- To obtain the cross 
section of LRG processes from the function R{Y — y2, z) we need to do four things: 



1. We have to insert the "hard" cross section of Eq. ( |2.30|) in all Pomeron lines at rapidity y2. 
That means that we have to change — > nz"~^ aj^t- 

2. To describe the absence of particle production in the interval Ay = yi — y2 we should 
replace 7 ^ 27 in the function R. As we have mentioned above, these rules are a direct 
consequence of the AGK cutting rules. Recall, that the factor Q in Eq. (|2.29|) can also be 



derived from the eikonal formula for the total cross section, making this substitution. 

3. To generate Pomerons at rapidity yi we should replace z in Eq. ( |2.33|) by 

z S{y2-yu2^,z) . (2.34) 

This replacement produces the generating function at the rapidity yi and the argument 27 
in Eq. ( p. 34]) , reflecting the fact that none of the Pomerons created in the interval from 2/2 
to yi, can be the source of the particle production. 
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4. Finally, we substitute S {yi, gpAih)) for z, since below rapidity yi we have normal "fan" 
diagrams for every Pomeron at rapidity level yi. 

The resulting formula is 

L{Y,yuy2,bt) = (2.35) 

' dR{y - y2,z,2-fy 
dz 

which can be rewritten in the form 



■ (^jet ■ S {yi,bt,'j) 



L{Y,y,,y2,h) = (2.36) 

9PN ■ 9PA{bt)e^^''-y'^y''^ {2KA{y2,bt) - KA{yubt) + lf 

{KA{yubt) + l) ■ {2KAiY,bt)-KAiyi,h) + l)' 

To find the survival probability of the LRG we need to divide Eq. ( p.29| ) by the inclusive cross 
section for di-jet production. We will show in the next subsection that this cross section is 

(Tinci [Y, yi, 1/2, bt) = a jet 7 — 7 , x , • (2.37) 

[K'AKyi.bt) + 1) 

Finally, we obtain for the survival probability of a LRG process the following formula 

< \S {Y,yi,y2) > = . „ , , „ , 777 , ,M • 2.38 

/ d^ b't d^ bt (Tinci {Y, yi, ?/2, h, b't) 

2.5 Inclusive production 

In this subsection we discuss the inclusive production of particles for hadron-nucleus interactions 
at high energy. It is well known that for inclusive production we have a number of very useful sum 
rules which follow from the AGK cutting rules |^| These sum rules simplify our calculations. 



reducing them to estimates of several Mueller diagrams [B^ . 



2.5.1 Single inclusive cross section 

The single inclusive cross section can be described by one Mueller diagram shown in Fig.|Tl] 



Recall, that the diagram of Fig.|TT] is what remains after the AGK cancellation in a system of 
Reggeon diagrams which includes "fan" diagrams in the entire region of rapidity {Y) and their 
eikonal repetition. The final result is 



Ay 



'^^^t±A = [ d'hgpMa^pe^^'^-y^SiyA) = ! d\ gp^ gpA{h) , (^.39) 
dy J J [KAiy, bt) + 1) 

where ap denotes the vertex of emission from the Pomeron. 
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Figure 11: The Mueller diagram for single inclusive cross section. 



2.5.2 Rapidity correlations 

For the double inclusive cross section, the only diagrams that survive the AGK cancellation are 



shown in Fig.|T2| and their contributions are 



dyi dy2 



2 

'9PN 



The rapidity correlation function is defined as 



1 



R{yi,y2) 



crtot dyi dy2 
1 dcr 1 d(j 
crtot dyi crtot dy2 



(2.40) 

(2.41) 
(2.42) 
(2.43) 



(2.44) 



Note that correlation function R is not equal to zero. In the case of the Glauber approach 
the second term in Eq. ( p.40|) ( see Eq. (|2.42|) ) cancels with 1 in Eq. (|2.44|) and only the first 
term ( Eq. ( p. 411 ) ) remains to generate a correlation function. In our approach we do not 
have such a cancellation and there are many correlations in the "fan" diagrams alone, without 
Glauber rescatterings. We present our numerical estimates for the correlation function in the 
next subsection. 



19 



yi 




Figure 12: The Mueller diagram for double inclusive cross section. 

2.6 Numerical estimates 

Proton - proton interaction: We first estimate phenomenologically the vertices of tlie Pomeron- 
Pomeron and tfie Pomeron-proton interactions using data from total, elastic and diffractive cross 
sections for high energy proton-proton collisions. Our goal is to describe proton-proton scattering 
using a phenomenological model as close as possible, to one Pomeron exchange which success- 
fully described the hadron data @]. However, we also want to include the experimental data on 
diffractive production which cannot be described in the one Pomeron exchange model. This is 
the reason why we decided to include " fan" Pomeron diagrams of Eq. ( p.l4|) in the proton-proton 
interaction. It is difficult to justify apriori such an approach, but developing such a model will 
allows us to extract the values of parameters for nucleon-nucleon interaction which we need to 
estimate the hadron-nucleus interaction. For the Pomeron-proton vertex we use Eq. ( |1.6| ) with 
A = 1 and we describe the the total cross section using the following generalization of Eq. ( |2.14] ): 

= (2.45) 

J J HN[y,bt) + 1 

We introduce similar modifications in Eq. ( |2.1| ) - Eq. ( |2.3| ), as well as in Eq. ( |2.27| ) and Eq. ( |2.28| ). 
Fitting the experimental data we obtain the following set of parameters: 

1. Rl = 25 GeV-'^ . 
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2. g%j,{bt = 0) = 70 GeV-^ 



3. 7 = O.UgpN{bt = 0) = 1.19 GeV^ . 



4. A 



0.07. 



Fig. shows our fit compared to the experimental data, where the survival probability for the 
case of p-p interaction was calculated using the eikonal approach. The results are reasonable for 
a first attempt. 
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Figure 13: Total ( Fig. 13-a ), elastic ( Fig.l3-b ), diffractive dissociation ( Fig. 13-c ) cross sections 
and the survival probability of LRG processes ( Fig. 13-d ) for proton - proton scattering. 

Proton - nucleus interaction: Using the parameters mentioned above, we calculate the total 
cross sections for the hadron- nucleus interaction at high energies, (see Fig.|l4| ). 



The corresponding elastic cross sections are shown in Fig. 15 



In Fig.|T6| we present the total diffractive production cross section as a function of energy. The 
single diffractive dissociation cross section, Eq. (|2.28|) , is presented in Fig.|l^ at fixed rapidity, 
y = In (S/Sq), ^/S = 2000 GeV and a/Sq = 1 GeV , where S denotes the center mass energy 
squared, as a function of the rapidity gap, y = ln(s/so), from to 15.2, which corresponds to 
= 1 - 2000 GeV and = 1 GeV. 



In Fig. |I8| we show our estimates for survival probabilities of large rapidity gaps for a jet produced 
in the interval of rapidity 10-12.5, as a function of energy. 
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Fig. 14-c Fig.l4-d 

Figure 14: Total proton - Nucleus cross sections for Ne ( Fig.l4-CL ), Mo ( Fig.l^-h ), Au ( Fig. 14-c ) 
and A=300( Fig.l^-d ). 



Estimates for single inclusive cross section, = dind/ cr total are given in Fig. |T9| as a function of 
rapidity, |/=0-15.2. The value of the correlation function in Fig. |! 
value of rapidity y2 = 10 as a function of yi from 10 to 15.2 for the fixed energy y/s 



is calculated for the fixed 
2000 GeV. 



3 Nucleus - Nucleus interaction 

3.1 Classification and selection rules for Reggeon diagrams 

Describing the nucleus-nucleus interaction is much more difficult than the description of hadron- 
nucleus scattering. In Fig. ^ one can see four topological classes of Reggeon diagrams, which are 
second order in Gsp, the triple Pomeron coupling. In general, they cannot be reduced to "fan" 
diagrams, as can be seen for the first diagram, which we shall call a "net" diagram. 

Using our main parameter k^, the rules of our selection are well defined as we have discussed in 
section 2.1. Complications arise as we have two parameters: and k^^. The general type of 
diagram that we select is the diagram which gives a contribution of the order of gp_Aigp-A2 ' ■ 
^Ai ■ '^^2- This means that we do not consider only diagrams with Pomeron loops. The whole 
structure of the diagrams looks rather complicated and, we hope, that summing all of them will 
yield a new example of high energy behaviour which is a natural generalization of the Glauber 
approach. Unfortunately, our new approach is quite different and more difficult than a simple 
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Figure 15: Elastic cross sections for Ne ( Fig. 15-a ), Mo ( Fig.l5-h ), Au ( Fig. 15-c ) and A=300( 
Fig.l5-d ) for proton - Nucleus scattering. 



summation of "fan" diagrams. As an instructive example, of how rich the system of the diagrams 



is, we show in Fig. ^ the "net" diagrams of order C'l^. 



3.2 An instructive example: the G|p - order Pomeron diagrams 

In order to demonstrate both the technique of calculation and the classification of the diagrams, 
using our general approach, we consider the "net" diagrams of order G\p, shown in Fig. ^ In 
Fig. Wl we have F|: 



1. Four diagrams with the same structure as the first diagram and their total contribution to 
the amplitude is 

Agl . gl . ■ 7^ ■ '\ - A^e^^ ^ e'^y'y'dy'j . (3.1) 

2. Two diagrams with the same structure as the second diagram, their total contribution to 
the amphtude is 



2gl ■ gl ■ e^^ ■ 7' ■ (e^" - e^^^y + A^e^^^ e^^' y' dy' + e^^Ay 



(3.2) 



"'^To abreviate we use the notation g\ — gp-Ai and 52 = .9P-A2 
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Fig. 16-a 
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Figure 16: Total diffractive production cross section as a function of energy for Ne ( Fig. 16-a ), Mo ( 
Fig.l6-h ), An ( Fig. 16-c ) and A=300( Fig.l6-d ) for proton - Nucleus scattering. 
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Fig. 17-c Fig.l7-d 

Figure 17: Single diffractive dissociation cross section at fixed rapidity, Y = ln{S/So), -v/S* = 2000 
GeV and \/5o = 1 GeV , as a function of rapidity gap, y = ln(s/so), from to 15.2 for 
= 1 - 2000 GeV and^=l GeV, for Ne ( Fig.l7-a ), Mo ( Fig.ll-h ), Au ( Fig. 17-c 
) and A=300( Fig.l7-d ) for proton - Nucleus scattering. 
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Figure 18: Survival probability for Ne ( Fig. 18-a ), Mo ( Fig.l8-b ), Au ( Fig. 18-c )and A 
Fig.l8-d ) for the fixed jet produced in the interval of rapidity 10-12.5 for proton - 
scattering. 
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Figure 19: Density of produced hadrons for Ne ( Fig. 19-a ), Mo ( Fig.l9-b ), Au ( Fig. 19-c ) and 
A=300( Fig.l9-d ) as a function of rapidity, Y, from to 15.2, for proton - Nucleus scat- 
tering. 
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Fig. 20-a 
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Figure 20: Correlation function for Ne ( Fig. 20-a ), Mo ( Fig.20-h ), Au ( Fig. 20-c ) and A=300( 
Fig.20-d ) for the fixed value of rapidity y2 = 10 as a function of yi from 10 to 15.2, for 
the fixed rapidity Y = ln{S/So), S = 2000 GeV and Sq = 1 GeV for proton - Nucleus 
scattering. 
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Figure 22: The "net" Pomeron diagrams for nucleus-nucleus intereaction of order G^, 
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3. Two diagrams with the same structure as the third diagram and their total contribution to 
the amphtude is 

2gl . gl . ■ 7^ ■ ^\ ^ + e^^Ayj . (3.3) 

4. Two diagrams with the same structure as the fourth diagram, their total contribution to 
the amplitude is 

2^7? ■ g'l ■ e^^ ■ 7' ■ (e^^ - e^""^ + A^e^y e'^y'y'dy' + e^^Ay) . (3.4) 

5. Two diagrams with the same structure as the fifth diagram, their total contribution to the 
amplitude is 



2g! . gl . e^y ■ 7^ ■ " ^ 



y] . (3.5) 



Summing all contributions given by Eq. ( |3.1| ) - Eq. ( |3.5| ) we obtain 

- gl ■ gl ■ e'^y ■ 7=^ ■ (3e=^^^ - lle^^^ + 9e^^ - 1 + 4Ae'^^i/) , (3.6) 



where the minus sign in front refiects the odd number of Pomeron loops in the diagrams. Eq. ( p.6| ) 
looks rather complicated. The key observation is that adding four digrams of "fan tree" type ( 
see the fourth diagram in Fig. 20 for the G|p- order and Fig. ^ for the general structure of these 
diagrams ) we obtain a very simple formula. Indeed, the four diagrams of the "fan tree" type 
give the contribution: 

- gl ■ gl ■ e'-y ■ 1' ■ (2e^^^ - 2 - 4e^^Ay) , (3.7) 



which leads to sum 

J2{Fig.20-l - Fig.20-5) + 4-Fig.20-A = 3 ■ gl ■ gl ■ e^y ■ -f^ ■ [e^y - if . (3.8) 

We, therefore, obtain the result in a very compact form. On the other hand the "fan tree" 
diagrams have a very simple form and their contribution can be calculated easily just by using 
the simple formula ( see Eq. ( ^.13| ) ) for the "fan" diagrams. 

To obtain the complete contribution of order Ggp we have to add to Eq. ( p.8| ) the contributions 



of the "fan" diagrams of the type given in Fig. 23 and subtract the contribution of Fig. 20-4 



diagrams. The contribution of the "fan" diagrams with "fan" looking down consists of two terms: 

1. One diagram of Fig. 21-1 which leads to 

lgi-gl-e^y-^'-{e^y-iy . (3.9) 

2. Four diagrams of Fig. 21-2, which give 

l^gi ■ g'2 ■ e'^y ■ 1' ■ {e'^y - if . (3.10) 
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The total contribution of these diagrams is 

-9i-9t-e''y-l'-{e^y-lf . (3.11) 

Summing Eq. ( p.8|) and Eq. ( p.ll|) we obtains a beautiful result 0: 

All contributions of G\p — order + 4 ■ Fig. 20 — 4 = (3.12) 

From Eq. ( |3.12| ) one can see that we can easily obtain the answer if we calculate the "fan tree" 
diagrams of Fig. |23 separately. 




Figure 24: The general structure of the "fan tree" diagrams. 

3.3 Sum of the Pomeron diagrams: general approach 

Based on the above example we have formulated a general algorithm: to the order of G^p the 
answer is 

i-ir 9i-92-e^y-^'-{e^y-iyig,+g,r - (3.13) 

n 

- ^ Ck{''fan tree" diagrams] C^p^Y) g\g2~^ , 

k=l 

^ We add the "fan" diagrams with the "fan" looking down to get Eq. ( |3.12|) . 
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where the contributions of the "fan tree" diagrams have to be evaluated. First, we discuss why we 
get such a simple first term in Eq. ( |3.13| ). In the simplest example of the diagrams in Fig. |21] one 
can see that if we integrate the "fan tree" diagram of Fig. ^T]-5 over the position of the upper G^p 
( ) we will get the contribution of the "net" diagram of Fig. pl|-4, for the region of integration 



yi < y2- However, this construction yields only one of the diagrams in Fig. ^T]-4. We can put it 
differently, that by integrating the diagrams in Fig. over yi and y2 without any restriction, 

we obtain a factor (^^{ ^ — 1 } ^ . i.e., we include two diagrams of Fig. 21-5 instead of one. 
Therefore, to obtain the correct answer we have to subtract one diagram of Fig. p]-4. 

Before formulating the rules as to what number of "fan tree" diagrams we need to subtract, let 
us consider in more detail the contribution of these diagrams (see Fig. ^ ). Replacing each line 
in the upper and low "fan" in Fig. by the full "fan" diagram function of Eq. (|2.13|) , we obtain 
the explicit expression for the sum of these diagrams: 



^{y,b,b') 



y g2 . ^2A(y-yi) . gA(j/i~j/2) 



(5-17 (e^(2^-fi) 



+ 



dyi 



91-92 2^ ^n+l,m+l (Gsp, 
n=l,m=l 



y)-9i 



yi 

'0 

92, 



9l ■ e^^^^ 



(^27(e^^^ + l) + l) 



Tdy2 



where 



9i 



9o,i 
TiRl 

Integrating over b and b' we have 



exp 



92 



90,2 

TiRi 



U2' 



m 



GIp ■ Rl ■ e^^ ■ vr^ . i?; 



/ e^^y-y^^dyi ■ / e^y^ ■ dy2 
Jo Jo 



gi 



X ■ dx 



'0 [x'j {e^(y-y^) 
The final expression for this diagram is 



i) + v. 



92 



y ■ dy 



GIp ■ RI- e 



■ TT^ . Rl 



dyi 



yi 



Ay 2 



dy2 



9il 



A(s/-s/i) _ ]^ 



^\7 {e^iy-y^) - 1) + 1 



Ln 



921 {e'^y 



921 [,e 



As/2 



Ln 



92l e 



1+1 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



where 



9i 



9o,i 
TiRi 



There are several important properties of the "fan tree" diagrams, which we will use below: 



1. These diagrams, in contrast to the "net" diagrams, yield a term, which is proportional to 
e^y. The coefficient in front of this term is 1 to all orders of G^p. 
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2. The first term in the expansion of these diagrams is proportional to the G|p ■ ■ 92- "Fan" 
diagrams also have a term e^^, but only have one of the vertices gi or g2 to the first power, 
for all orders of G^p. 

3. As has been stated above, to obtain a compact expression for the sum of the diagrams in 
the nucleus-nucleus amplitude at high energy to all orders of G^p we need to add a definite 
number of "fan tree" diagrams of the same order. 

We now begin our calculation of the amplitude using these properties of the "fan tree" diagrams. 

First of all, we introduce a function, which contains all the "fan" diagrams, all the "net" diagrams 
and part of the "fan tree" diagrams, namely, those contained in the first term in Eq. ( p.l3 ): 

the term of order 7^ in Eq. (|3.18| ) is Eq. ( p.l2|) . To find the second term in Eq. (|3.13|) we expand 
Eq. ( p.l8| ) with respect to the powers of 7: 



1. To first order in 7 we have: 

-g,.g,.e^y-^-ig,+g2){e^y-l) . (3.19) 



This is the first contribution from the "fan" diagrams. Fig. ^ 
2. Second order in 7 has the form 

2■g^■g2■e^y■^'■ (g, + ^72)' (e^^ - l)' • (3.20) 

One can see that this term which is proportional to 2'y'^gf ■ g^ ■ e^^ in Eq. ( |3.20| ) has the 
same structure as a term of the same order in the "fan tree" diagram in the total amplitude. 
The difference is only in the value of the numerical coefficient: in the amplitude this term 
has coefficient 1 instead of the 2 in Eq. ( p.20| ). Therefore, using the properties of the "tree" 
diagrams we can conclude, that to obtain the exact coefficient 1 we must subtract from this 
term the first term of the expansion Eq. ( p.l4| ), namely, g1- g^' ^2,2 (Gsp, y), which also has 
a coefficient 1. 



3. The third order in 7 is 



g,.g^.e^y.^^.(e''y-l)\g, + g2f. (3.21) 



In order to obtain a coefficient 1, we subtract from this term the second term of the expansion 
Eq. ( CT) multiplied by 2, 



2-[gl-gl-A2,^ + gl-gl-A^^2) ■ (3.22) 
This gives us four diagrams of the required order. 
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Conforming this expansion we find that the number of "fan tree" diagrams of order ■ ■'j"''^"^ ^ 
which we have to subtract is equal to 

(n + m — 2)! 



(n- 1)! (m - 1)1 
This means that the function to be subtracted is equal to 



1 . (3.23) 



F2{y,b,b')=g,-g2 E f ^^^^^ - iV ^n+i,™+i ■ < ■ ■ (3-24) 
Using Eq. ( p.l4| ), we obtain 



F, {y, b, b') = 2- Glp ■ gl ■ gl ■ e^^ • / e^^^-y^^dy^ ■ / e^^^ ■ dy^ ■ 

Jo Jo 

^ o -$(?/, 6, 6') . (3.25) 

{g,^(e^iy~yr)-l)+g,^(e^y2-i) + if 

Finally, the two particle irreducible set of Pomeron diagrams which gives us the opacity Q{s,bt) 
( see Eq. ( |TT^ ) and Eq. ( |rT3|) ) is equal to 

n{s,bt) = J d% F{y,b,\bt~b\) , (3.26) 

where 

Fiy,b,b') = F,{y,b,b') - F^iy^h') = 

= $1(1/, 6, 6') + $2 (?/,&,&') - ^z{yAh') , (3.27) 

^,{y,b,b') = F,{y,b,b')= ' (^.28) 

(^1+^2)7(6^^ - 1) + 1 

$2(1/, 6, 6') = $(z/,6,&') (3.29) 

^ '""io ((7i7 {e^^yy^'^ - 1) + 1)' ^''^0 (^27 {e^y^ + i) + if 

$3 (y, &, b') = 2- GIp ■ gl ■ gl ■ e^^ ■ / e^^^-^^^t/yi ■ / e^^^ ■ dy, ■ (3.30) 

JO JO 



{gii {e^^y-y^) - i) + ^?27 (e^^^ - 1) + 1) 



3.4 The Ultra high energy asymptote 



3 • 



The asymptotic form at ultra high energy of the opacity Vt originates in Eq. (|3.27|) from the 



difference between the second and the third terms in this equation. To our surprise this asymptotic 
behaviour is just the exchange of the Pomeron but with an intercept which is two times smaller 
than the intercept of the input Pomeron: 

f(y.M'U„, = ||-^-e*-». (3.31) 
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The opacity Q{s, ht) is equal to 



^asymp 



[sA) = I d%F,„{yAh') = ^-^-e^y ■ ^! ' ■ exp"^^ . (3.32) 
^ ^ ^ Af + Al 



Using the general formulae of Eq. (|2.1|) - Eq. (|2.5| ) and Eq. ( p.32|) we can calculate the asymptotic 



behaviour of the total and elastic cross sections in the kinematic region where the opacity f2 is 
small ( < 1 ). These calculations lead to 

oTor\y) = 2 I d%tn^„{y,bt) Af x 4 x (3.33) 



5 



<r'iy) = I d\nly^,{yM) — > ^! ""-^i xe^^. (3.34) 
^ + Al 

We will discuss below how to calculate processes of diffraction dissociation in our approach, but for 
the sake of completeness we present here the result of the asymptotic behaviour of the diffractive 
cross section at ultra high energy. It turns out that this cross section is 

oc ^ X (A, + A2) X e^y . (3.35) 



It should be stressed that Eq. ( p.35|) is quite different from what we expect if the asymptote is 



a real Pomeron exchange with a different intercept. Let us recall that in the single Pomeron 
exchange model 



SD asymp 2Apy 

^asymp ^ el ^ ^ 



The formula for the inclusive cross section in the kinematic region of sufficiently small VL is also 



quite different from the single Pomeron exchange (see Fig. (|25|) ). Indeed, from Fig. (|25|) one can 
see that the single inclusive cross section results from the sum of "fan" ( or better to say "tree 
fan" ) diagrams. Therefore, the density of produced particles in a unit of rapidity 



I asymp i j 

^ ^ d<^ind /dye 

P — asymp 
^tot 



has an asymptotic limit 



1 1 



p oc A^~^ ■ A2' ■ e-'^'-y — y 0. (3.36) 



Therefore, we can conclude, that the Pomeron interaction with both nuclei leads to a small 
number of particles produced in the central rapidity region. This is the most striking difference 
of the resulting effective Pomeron from the input Pomeron, which has a uniform distribution of 
particles in rapidity. 

It is also interesting to note, that in our amplitude the second term has a stronger energy de- 
pendence, in spite of the fact, that each diagram in this term has a weaker energy dependence, 
than each diagram of the same order in the first term. This means, that the sum of the diagrams 
with the weak energy dependence grows with energy faster, than the sum of the diagrams with a 
stronger energy dependence. 
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Figure 25: The Mueller diagram for the inclusive cross section in our approach for a nucleus-nucleus 
collision. 
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3.5 What energy is asymptotic? 



As has been mentioned, the three terms in Eq. ( |3.27| ) have different asymptotic behaviour. At 
ultra high energy the last two terms dominate and they provide the leading asymptotic behaviour 
that we have discussed in the previous section. To estimate the value of energy that we can 
consider to be ultra high, so that we can safely apply our asymptotic formulae of Eq. ( |3.31| ) - 
Eq. ( p. 361) , we calculate the three terms of Eq. (|3.27|) separately, namely. 



ai = 2 y d%(fh' ■ $1 {y, b, b') 
£72 = 2 y d%d%' ■ $2 {y, b, b') 
(73 = 2 1 $3 (?/,&,&') 



where the total cross section at small Vl is equal to 



(^tot = 0-1 + 0-2-0-3 



(3.37) 
(3.38) 
(3.39) 

(3.40) 
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Figure 26: The Energy dependence of the first ( Fig.26-a ),the second ( Fig.26-h ) and the third ( 
Fig.26-c ) terms of Eq. \3.4-(\ ) as well as the value of the total cross section ( Fig.26-d ). 
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The result of the calculation is given in Fig. (|2^) where the parameters of section 2.6 were used 
for the proton-proton interaction. One can see from Fig.(p6D that in a range of energy, up to the 
Tevatron energies, the first term is much larger than the other two. Therefore, we can conclude: 



• That ultra high energies, where our asymptotic answer is valid, are above of the Tevatron 
energies, and we will only be able to see the characteristic features of the real asymptote at 
the LHC. 

• We can take into account only the first term for all numerical estimates for the RHIC 
energies as well as for all energies below the Tevatron energy. This is the reason for only 
taking the first term in the following numerical estimates. 



3.6 Total and elastic cross sections 



To calculate total and elastic cross sections we have to use the general formulae of Eq. and 
Eq. ( |2.2D with the opacity Q{s, ht) defined in Eq. ( |3.32D through the function 6, h') given in 
Eq. (p. 271). The results of these calculation are shown in Fig.(p7|). 



3.7 Diffractive dissociation processes 

We start the discussion of the diffractive cross section for nucleus-nucleus interaction by rewriting 
the formula for the amplitude in terms of the function S{y^ ht) which describes the sum of the 
"fan" diagrams ( see Eq. ( p.l3| ) ). We introduce two such functions: and which sum the 
"fan" diagrams with up and down looking "fans", respectively. Using these functions we rewrite 
our amplitude as 

n=QO 

F {y, 6, b') =g,.g,.e^y-S^-SrT.i^- ^iT (1 - ^t)" • (3-41) 

n=0 

This is a very convenient representation of the amplitude, because, we can use the AGK cutting 



rules and reduce the problem of finding the diffractive cross section to the problem of the 
diffractive cross section for hadron-nucleus interaction that has been solved in section 2.3. To see 
this, we rewrite the amplitude in the form: 

F (y, b, b') = ^1 . ^2 ■ e^^ ■ 5f E (1 - ^t)" E ("l)'" • (3-42) 

n=0 m=0 

This sum can be rewritten in terms of the amplitude for single diffraction dissociation. We 
calculate the diffraction dissociation of the nucleus A2. In this case we notice that the cut of Si 
leads to a multiparticle production process or to the diffractive dissociation of the nucleus Ai and 
, therefore, it does not contribute to the process of interest. Finally, we obtain 
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Figure 27: Total and elastic Ai — A2 cross sections versus energy. 
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n=0 m=0 



A:=l 



The function S^^^ differs from D{y,yi,bt) of Eq. ( |2.28|) only by some factors, namely, 
After simple algebra, we reduce Eq. (|3.44|) to 

n=oo m=n 



n=0 m,=0 



m+l 



fc=i k\{m-k + l) 

91 {b - b') ■ g, m ■ e^^ ■ {2S^ - ^) 
l-(l-5t)(l + 25i-Sf^) 

+ 1 - (1 - 5,) (1 - 25,) • ^'-^'^ 

Using Eq. ( |3.26| ) and Eq. ( |3.45| ) , we can obtain a simple expression for the cross section of the 
single diffractive production: 

a'^ = J d\ ■ (Ff ^ {y, h) + Ff ^ (y, h)) e'^^^'''') . (3.46) 

We can also obtain a simple formula for the total cross section of diffractive dissociation directly 
using the unitarity constraint of Eq. ( p. 221 ) which leads to 

F^^{yAb') = 2 9i-92-e^y 



91 + 92) 7 (e^J' - 1) + 1 



2 91-92-e'^y 

2{9^ + 92)i{e^y ~l) + 1 ■ ^ • ^ 



Eq. ( p.47[) leads to the total cross section of diffraction dissociation: 



j d\ F^^ e-^(J^'^*) . (3.48) 



The calculations were performed for Ai — A2 interaction with {Ai = 20) - {A2 = 100), {Ai = 
200) - {A2 = 100), {Ai = 200) - {A2 = 300). The results for the total diffractive dissociation 



cross sections are shown in Fig.^ The single diffractive dissociation in Fig.^ is calculated 
as a function of rapidity gap y for a fixed rapidity Y = \nS/ Sq for the ^/S = 2000 GeV and 
^/So = IGeV. All parameters were taken to be the same as for the elastic and total cross sections 
calculations. 
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12.5 13.8 14.6 15.2 

Figure 29: The single diffractive dissociaton Ai — A2 cross sections as a function of rapidity gap y for 
a fixed rapidity Y = hiS/So for the = 2000 GeV and ^/S^ =1 GeV. 



43 



3.8 Survival probability of Large Rapidity Gaps 

We can calculate the survival probability of the LRG processes using the expression for the 
scattering amplitude in the form of Eq. ( |3.41 ) which we rewrite as 

n=oo m=n k=n 

F {y, b, h') = g.-g^-e^y-Y^HCl {-S,r^' E (-^t)'^' • (3-49) 

n=0 m=0 k=0 

Substituting the function S^^ instead of one of in Eq. ( p.49| ) we obtain the survival probability 
for a process of di-jet production which can be accompanied by hadrons with rapidities which 
are smaller than 1/2. The function S^^ differs from L{y,yi,y2,bt), found in Eq. ( p.37| ), by some 
factors, 

^f^ = (3.50) 



2 



g-A(j/2-i/i) (^1 ■ 7 ■ (2e^y^ - e^y^ - l) + 1 

^i-* ((71 ■ 7 ■ {e^y^ - 1) + 1) ■ ((71 ■ 7 ■ {2e^y - e^v^ - 1) + 1)' ' 

Replacing one of the by Sf^ and performing the explicit summation over all the indices, we 
obtain: 



(,,(,') = g.-m-e'^'. flf'^'\ (3,51) 

^ . . . Sr ■ + ,3,2) 

' ((91 + 92) 7 (e^» - 1 ) + 1)' ' 

The total cross section for the two particle irreducible diagrams is a sum of two terms F^^ = 
F^^ + F^^ which is equal to 



^''{y,b,b') = g^ ■ 92 ■ ■ {Sf^ + Sf^) 

(ga (e^y - 1) + 1) (921 (e^y 



(3.53) 



((^7i+^72)7(e^^-l) + l)' 
where Sf^ is given by Eq. ( p.50| ) and 

^f^ = (3.54) 

^-A{y2-yi) (92- 1- {2e^^y-y'^ - e^^y-y^^ -i^ + i^ 



{92-1 ■ [e^^y-y^^ - 1) + 1) {92- 1- {2e^y - e^(y-y^) - 1) + 1)^ 



For the calculation of the survival probability we also need the expression for the inclusive process 
in the case of nucleus-nucleus interactions. 
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Figure 30: The Mueller diagram for the inclusive process for a A-A interaction. 



The Mueller diagram for the inclusive cross section is given in Fig.^ and the analytic expression 
is 

Finally, from Eq. (|2.32| ) we obtain 

< " > = ' ^'-''^ 

where f2 is defined by Eq. ( ^.261 ). Using Eq. ( p.5(j| ) we performed calculations with [Ai = 20) 
- {A2 = 100), {Ai = 200) - {A2 = 100), (Ai = 200) - {A2 = 300), for di-jet production in the 
rapidity interval 10-12.5. The result of the calculation is given in Fig.^ as a survival probability 
versus energy. 



One can see in Fig. 31 that the value of the survival probability turns out to be rather small. 



3.9 Inclusive production 

In this subsection we repeat the calculation, that has been performed for hadron-nucleus collision 
in section 2.5. 
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Figure 32: The multiplicity of produced hadrons for Ne-Mo, Mo-Au, Au-A=300 in Ai — A2 interaction. 
3.9.1 Single inclusive cross section 

The Mueller diagram is shown in Fig. ^ This diagram leads to a simple formula 

da{Ai + A2) 



dye 



J cfb d%' ap Sa^ {y - yc, b) ■ (yc, b') 



(3.57) 



( f^A^ iv-Vcb) + 1 ) ( {yc, b) + 1) 



Integrating this equation over yc and dividing by the atot we obtain the multiplicity N of the 
produced hadrons in A-A interaction, where the results are shown in Fig.|32| for the cases of 
interactions of Ne-Mo, Mo-Au, Au-A=300. 

It is interesting to note that the density of produced hadrons in the central rapidity region at 
high energy has a very simple relation to the density in nucleon-nucleon collision: 



da{Ai+A2) 
dye 

atot{Ai + A2] 



g'p_j, AI ■ InAl-lnAI 



7 



1 + (*) 



da{N+N) 
djjc 

atotiN + N) 



(3.58) 
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where we assumed 

R\^R%- (3.59) 

and 

A2 > Ai. (3.60) 

3.9.2 Rapidity correlations 

Fig.^ shows the Mueller diagrams for the double inclusive cross section which can be written 
analytically as 



^l^(^l±M = ^a^f [ d% d%' . (3.61) 
dyi dy2 J 

.^Hy2~y.) (^ _ h) Sa, (2/1, h') (3.62) 

+ Sa^ (y - y2, b) Sa^ (2/2, b') Sa2 (y - yi, b) j d%" Sa2 (yi, b") (3.63) 

+ Cii/'e^(2,'-,,-,,) _ ^) ^y^^ ^^^^ ^^ ^^^^ 

+ G,p / dy'e'^^y^+y^-^y ^Sa, (y', b') Sa, {y - y2. b) Sa, {y-yi,b) , (3.65) 
Jo 



where each term of Eq. ( |3.62| ) - Eq. ( p.65| ) corresponds to Fig. - Fig. ^4. The asymptotic 
behaviour for the double inclusive cross section is: 



dMA + A^) (a^)2 



dyi dy2 7 



2 



{rI,R%) \n^4 In' 4 . (3.66) 



4 Conclusions 

In this paper we found a natural generalization of the Glauber approach to hadron-nucleus and 
nucleus-nucleus collisions. It has been known for three decades |1^ that such a generalization 
requires finding a theoretical way to take the triple Pomeron interaction into account or, in 
other words, to include the diffractive dissociation processes. In the framework of the Pomeron 
approach, this problem was solved here using two main physical ideas, or more precisely two small 
phenomenological parameters: 

1. Typical distances for the Pomeron structure which are sufficiently small, justify using the 
pQCD approach for the evaluation of different types of the Pomeron interaction. 

2. A small value of the triple Pomeron vertex {G^p/ qp^n ^ 1) which allows us to use the 
following set of parameters in summing the Pomeron diagrams: 

G3P . 1 

7 = — — ^ 1 , 7 X As ^ 1 . 
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Figure 33: The Mueller diagram for the double inclusive process in an Ai — A2 interaction. 

Based on these parameters we were able to formulate the selection rules for the Pomeron diagrams 
and sum them. This lead to a significant generalization of the oversimplified eikonal type model 
for shadowing corrections without losing the advantages of such models. In this approach we 
calculated diffractive dissociation processes and survival probability of the large rapidity gap 
processes, and calculated the shadowing (screening) corrections for a large class of the "soft" 
interaction processes at high energies. 

The physics of our generahzation is very simple. In the parton model, the Glauber (eikonal) 
approach takes into account only the interaction of the fastest parton with the target. In our 
approach, we consider the interactions of all partons both with the target and with the projectile, 
if the last is a heavy nucleus. We are successful in finding a simple closed expression for such 
types of interaction making use of the fact that the nucleus is a dense parton system. 

Our main results are: 

1. We have generalized the Glauber approach for the nucleus-nucleus interaction at high en- 
ergies. 

2. We show that the asymptote of the two particle irreducible diagrams leads to the exchange 
of an effective Pomeron with an intercept which is two times smaller than the intercept of 
the input Pomeron. 
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3. We show that the asymptotic behaviour for nucleus-nucleus collision starts at sufficiently 
high energies ( higher than the Tevatron ). 

4. A systematic approach to diffractive dissociation has been developed both for hadron- 
nucleus and nucleus-nucleus collisions. 

5. Our approach leads to saturation of the rapidity density of produced hadrons 



MA^+^^) ^2 4t.in4^.1n/l^ ^^^^^ 

dye I yP-N ^2 f-i-^l in/l2 dye 



(4.67) 



This prediction is quite different from the Glauber-Gribov result given in Eq. ( p..l7| ). In 
particular, Eq. (f4.67|) predicts that the multiplicity for nucleus - nucleus collisions < N > 



{Ai - A2) = A't < N > {N ~ N) where Ai < A2, while Eq. ( |1.17|) leads to < iV > 

2 

{Ai — A2) = Ai ■ A2 . We predict the same dependence also for central collisions. 

On the other hand, the prediction of Eq. ( |4.67| ) is quite different from the predictions of 



the high density QCD approach u^. Indeed, the high density QCD approach leads to a 



formula which is very similar to Eq. ( 4.67 ), but due to the integration over the transverse 

momentum of the emitted partons the extra factor appears in front of the equation. 

2 

In particular, it turns out that < N > {Ai - A2) oc Af ■ Q'j^{s,A)- < N > {N - N). 
Q'^{s,A) is the saturation scale which increases with A and energy squared s. We expect 



that Qlis, A) grows with A ( at least Ql oc Ai [|T2| but it could even be proportional to As 



2/'o A\ ^T"^T,ro T,ri^-V, A f f^i- /n2 ^ ||x uuu III uuulu even uc piupui Liuiiai 

2 

38[| ), and , therefore, high density QCD predicts < N > {Ai — A2) oc A^ x 



2 



where A2 > Ai. 



^2 



6. The survival probability of the LRG processes has been calculated both for the hadron- 
nucleus and the nucleus-nucleus collision. 

7. The theoretical solution to the nucleus-nucleus interaction can serve as guiding light for a 
theoretical approach to the dense parton system. 

8. The numerical estimates which we have made, serve as a basis for run of the mill physics 
estimates. Hence, results found at the next generation of accelerators such as RHIC and 
LHC, should be compared to our estimates to determine if there are any unusual features. 

The most important corrections to our approach stem from the enhanced diagram ( see Fig.^ ). 
Direct calculation of the ratioQ 

Enhanced diagram G^p 2 



^^Fan" diagram 2 a'p 



r (4.68) 



^We thank K. Tuchin for help m calculating Eq. ( |4.68| ). 
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gP-A A 

"Fan" diagrams First enhanced diagram 

Figure 34: The first correction to our approach. 

shows that in our approach we have achieved reasonable accuracy only for sufficiently small values 
of the triple Pomeron vertex. In pQCD (see Eq. (|1.8| ) and Eq. ( |1.9|) ) one can see that the ratio of 
Eq. ( [4.68| ) is proportional to and, therefore, it appears that we can neglect the contribution 

of the enhanced diagram in the wide region of energy y = Ins < 1/as- However, the numerical 
parameters of our fit (see section 2.6 ) gives Glp/2a'p ~ 1. 
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